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I .  INTRODUCTION 


Neutron  damage  in  silicon  plays  a central  role  in  radiation  effects 
of  military  interest,  and  a knowledge  of  the  energy  dependence  of  this 
neutron  damage  is  required  to  relate  threat  and  simulator  spectra  as 
well  as  to  provide  a basis  for  a 1-MeV  equivalence  definition  for  hard- 
ness assurance.  To  meet  this  requirement,  a model  has  been  formulated 
to  effectively  utilize  all  available  neutron  cross  section  information 

1 2 

to  calculate  the  energy  available  for  displacement  damage,  ’ and  the 
calculated  results  have  been  compared  with  experimentally  determined 

3 

damage  (lifetime  degradation)  . In  this  report  calculated  and  experi- 
mentally measured  results  for  both  threat  and  simulator  spectra  will 
be  presented. 

Neutron  irradiation  of  silicon  produces  changes  in  the  electrical 
properties  of  the  material,  and  these  changes  manifest  themselves  in 
two  ways,  as  a degradation  of  minority-carrier  lifetime  and  as  a reduc- 
tion of  the  equilibrium  charge-carrier  concentration.  Radiation-induced 
neutron  damage  causes  permanent  changes  in  the  operating  characteristics 
of  silicon  devices  and  the  extent  of  this  permanent  damage  is  neutron- 
energy  dependent.  The  evaluation  of  energy  dependence  is  complicated 
by  the  annealing  of  unstable  damage,  by  effects  peculiar  to  devices 
if  devices  are  used  in  the  evaluation,  and  by  the  superposition  of 
lifetime  and  carrier-concentration  effects  if  both  change  significantly. 
Past  research  by  many  investigators  has  produced  substantial  information 
on  energy  dependence  and  this  will  be  discussed  in  the  following  para- 
graphs of  this  introduction.  Then  an  experimental  technique  for  measur- 
ing lifetime  degradation  without  interference  from  carrier-removal  and 
a calculational  procedure  that  utilizes  all  available  neutron  cross 
section  information  to  determine  the  energy  available  for  displacement 
damage  will  be  described,  and  both  calculated  and  measured  results 
will  be  presented. 

Calculation  of  displacement  damage  as  a function  of  neutron  energy, 
with  consideration  of  the  different  reactions  and  angular  distributions 


1.  J.E.  Youngblood,  W.R.  Van  Antwerp,  and  R.M.  Tapphorn,  "Displacement 
Damage  in  Silicon  Irradiated  with  6-  to  10-MeV  Neutrons",  USABRL 
Memorandum  Report  No.  2738,  April  1977.  (AD  #A039774) 

2.  J.E.  Youngblood  and  W.R.  Van  Antwerp,  "Calculated  Energy  Dependence 
of  Neutron- Induced  Displacement  Damage  in  Silicon",  USABRL 
Memorandum  Report  No.  2759,  June  1977.  (AD  #041465) 

3.  J.E.  Youngblood,  C.E.  Hollandsworth,  and  W.R.  Van  Antwerp,  "Neutron 
Damage  in  Silicon  From  Neutrons  With  Energy  Near  1-MeV",  USABRL 
Memorandum  Report  No.  2768,  July  1977.  (AD  #A043268) 
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of  the  neutrons,  evolved  through  the  work  of  Smith  et  al.,  Stein, 

and  Holmes.  These  efforts  and  the  more  recent  work  of  Guenzer  and 
7 8 

Manning  and  of  Rogers  et  al.  determine  the  fraction  of  energy  availa- 
ble for  displacement  damage  using  the  stopping-power  theory  of  Lindhard, 

q 

Scharff,  and  Schiott  (LSS).  The  LSS  theory  provided  a more  realistic 
basis  for  energy  partition  than  earlier  theories  and  it  was  experi- 
mentally confirmed  (by  measuring  the  ionization  fraction)  by  Satler 

and  Vook.  In  efforts  to  proceed  beyond  the  energy  available. 

Holmes,6  Curtis,1*  and  Gregory1 ^ have  attempted  to  analyze  in  detail 
the  electrical  behavior  at  disordered  regions.  Both  point  and  cluster 
defects  have  been  considered  and  a recent  review  of  work  in  this  area 

13 

has  been  presented  by  van  Lint  and  Leadon.  Cluster  formation  and 


4.  E.C.  Smith,  et  al.,  "Theoretical  and  Experimental  Determinations 
of  Neutron  Energy  Deposition  in  Silicon",  IEEE  Trans.  Nucl.  Sci., 

NS- 13,  No.  6,  (1966).  D.  Binder,  et  al.,  "Analytical  and  Experi- 
mental Predictions  of  Fusion  Neutron  Radiation  Effects",  AFWL-TR- 
66-41,  Vol.  1,  Hughes -Fullerton,  (1966). 

5.  H.J.  Stein,  "Energy  Dependence  of  Neutron  Damage  in  Silicon", 

J.  Appl . Phys.,  38,  No.  1,  (1967). 

6.  R.R.  Holmes,  "Energy  Dependence  for  Carrier  Removal  and  Lifetime 
Damage  by  Fast  Neutrons  in  Silicon",  Bell  Telephone  Laboratories 
Weapons  Effects  Studies,  Report  to  ABMDA,  Vol.  II,  Supplement  III, 
pp  67-88,  October  1,  1970. 

7.  C.S.  Guenzer  and  Irwin  Manning,  "Calculated  Displacement  Damage  by 
Neutrons  in  InSb",  IEEE  Trans.  Nucl.  Sci.,  NS-21,  No.  6,  (1971). 

8.  V.C.  Rogers,  L.  Harris,  Jr.,  D.K.  Steinman,  and  D.E.  Bryan,  "Silicon 
Ionization  and  Displacement  Kerma  for  Neutrons  from  Thermal  to 

20  MeV",  IEEE  Trans.  Nucl.  Sci.,  NS-22,  No.  6,  (1975).  Also, 

Op.  Cit.,  NS-23,  No.  1,  875,  (1976). 

9.  J.  Lindhard,  M.  Scharff,  and  H.E.  Schiott,  "Range  Concepts  and 
Heavy  Ion  Ranges",  Kg.  Danske  Videnskeb  Selskab,  Mat.  Fys.  Medd. , 

33,  No.  14,  1-42,  (1963). 

10.  A. R.  Sattler  and  F.L.  Vook,  "Partition  of  the  Average  Energy 
Deposited  in  Silicon  as  a Function  of  Incident  Neutron  Energy", 

Phys.  Rev.,  L55,  No.  2,  211  (1967). 

11.  O.L.  Curtis,  Jr.,  "Statistics  of  Carrier  Recombination  at  Disordered 
Regions  in  Semiconductors",  J.  Appl.  Phys.,  3£,  3109,  (1968). 

12.  B.L.  Gregory,  "Minority  Carrier  Recombination  in  Neutron- Irradiated 
Silicon",  IEEE  Trans.  Nucl.  Sci.,  NS-16,  No.  6,  53-62,  December 
(1969). 

13.  V. A.  J.  van  Lint  and  R.E.  Leadon,  "Implications  of  Cluster  Model 
of  Neutron  Effects  in  Silicon",  Lattice  Defects  in  Semiconductors, 
1974,  Inst.  Phys.  Conf.  Ser.,  No.  23,  The  Institute  of  Physics, 
London  and  Bristol,  pp  227-232,  (1975). 
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annealing  is  receiving  considerable  current  interest  because  of  its 

relation  to  fusion-energy  studies.  Recent  results1*1'1^  imply  that  there 
is  a threshold  energy  for  the  formation  of  stable  clusters,  that  there 
is  much  diffusion  out  of  clusters,  and  (thus)  that  there  is  a limited 
significance  to  clustering.  In  the  calculations  on  which  the  present 
results  are  based,  the  damage  is  modeled  only  to  the  extent  of  deter- 
mining the  fraction  of  the  energy  available  for  displacements,  and 
emphasis  has  been  placed  on  the  detailed  use  of  all  available  neutron 
cross  section  data. 

Experimental  work  on  the  energy  dependence  of  neutron  damage  has 
been  very  limited.  Work  with  monoenergetic  neutrons  has  been  done  by 

Smits  and  Stein,16  Cleland  et  al.,17  and  Speers.1®  Working  with  selected- 

19 

average  energies  from  polychromatic  sources,  Kantz  and  Lohkamp  and 
20 

McKenzie  have  obtained  experimental  results  related  to  energy  dependence. 

Smits  and  Stein16  compared  the  damage  in  silicon  produced  by  14  MeV 
neutrons  relative  to  reactor  neutrons.  They  found  a ratio  of  3.0  for 
carrier  removal  and  a ratio  of  3.5  for  charge  carrier  lifetime  degrada- 

2i 

tion.  One  should  note,  however,  that  Coppage"  has  shown  that  errors  of 


14". Yu.  V.  Martynenko,  "Annealing  and  Clustering  of  Defects  in  Cas- 
cades", Rad.  Effects,  29,  192,  (1976). 

15.  V.L.  Vinetskii  and  A.V.  Kondrachuk,  "The  Kinetics  of  Formation  and 
the  Parameters  of  Radiation  Defect  Clusters  in  Silicon",  Rad. 
Effects,  30,  227,  (1976). 

16.  F.M.  Smits  and  H.J.  Stein,  "Energy  Dependence  of  Neutron  Damage  in 
Silicon-Experimental",  Bull.  Am.  Phys.  Soc.  9,  No.  3,  238,  (19641. 
F.M.  Smits,  "On  the  Energy  Dependence  of  Neutron  Damage  in 
Semiconductors",  Sandia  Report  No.  SC-R-64-196,  (1964). 

17.  J.W.  Cleland,  R.E.  Bass,  and  J.H.  Crawford,  Jr.,  "The  Nature  and 
Yield  of  Neutron- Induced  Defects  in  Semiconductors",  Conference 

on  Radiation  Damage  in  Semiconductors,  Paris,  (1964),  Proc.  of  the 
7th  Int.  Conf.  on  the  Physics  of  Semiconductors,  Vol.  3,  Radiation 
Damage  in  Semiconductors,  t’aris-Royaument  1964,  pp  401-406, 
Academic  Press,  New  York,  (1965). 

18.  R.R.  Speers,  "Neutron  Energy  Dependence  of  Excess  Charge  Carrier 
Lifetime  Degradation  in  Silicon",  IEEE  Trans.  Nucl.  Sci.,  NS-15, 
No.  5,  9-17,  (1968). 

19.  A. D.  Kantz,  "Average  Neutron  Energy  of  Reactor  Spectra  and  Its 
Influence  on  Displacement",  J.  Appl.  Phys.  3£,  No.  7,  1944-1952, 
(1963). 

20.  J.E.  Lohkamp  and  J.M  McKenzie,  "Measurement  of  the  Energy  Depend- 
ence of  Neutron  Damage  in  Silicon  Devices",  IEEE  Trans.  Nucl.  Sci. 
NS- 22,  No.  6,  2319-2325,  December  (1975). 

21.  N.  Coppage,  "The  Influence  of  Dosimetry  on  Earlier  Damage  Equi- 
valence Ratios",  IEEE  Trans.  Nucl.  Sci.,  NS-22,  No.  6,  2336-2339, 
December  (1975). 


30  percent  are  likely  in  early  fluence  measurements  at  reactors  and 
these  would  result  in  equal  errors  in  quoted  damage  ratios.  Cleland 

et  al.^7  used  monoenergetic  neutrons  from  an  accelerator  but  based  their 
conclusions  on  observations  of  carrier  removal  and  thus  required  high 
fluences  and  attendant  long  runs.  The  present  experimental  work  is 

IS 

most  closely  related  to  that  of  Speers,  but  uses  different  readout 
techniques  and  better  fluence  measurements. 


II.  EXPERIMENTAL  TECHNIQUE 

The  experimental  measurements  of  damage  were  done  with  wide- base 
silicon  diodes  operating  at  a fixed  forward  current  of  0.1A  to  maintain 
a constant,  high  level  of  charge  injection  and  a high  sensitivity  to 
neutron  damage.  All  diodes  were  made  from  100  ohm-cm  p-type  silicon 

2 

and  had  a cross-sectional  area  of  about  0.1  cm  . However,  two  base 
widths,  about  1.27  mm  (0.050  inch,  "50-mil")  and  about  0.76  mm  (0.030 
inch,  "30-mil"),  were  used.  In  both  cases  a constant  injection  level 
existed  because  the  low  level  of  neutron  exposures  used  changed  the 
injected-carrier  lifetime  without  significantly  changing  the  equilibrium 
carrier  density.  Observations  were  made  on  the  change  in  forward 
voltage  (with  fixed  current)  induced  by  neutron  irradiations.  We  have 

shown  this  is  a direct  evaluation  of  lifetime  degradation.  The  50-mil 
diodes,  and  the  readout  technique  used  with  them  here,  were  used  exten- 
12  3 

sively  in  past  work  ’ ’ and  this  use  will  be  described  next.  Special 
efforts  were  made  with  the  30-mil  diodes  to  control  annealing  and 
injection  effects  and  these  will  be  described  at  the  end  of  this  section. 

An  analysis  of  the  current- voltage  characteristics  of  wide-based 
silicon  diodes  and  the  effects  of  neutron  irradiation  can  be  found  in 

22 

the  paper  by  Swartz  and  Thurston.  “ They  combine  the  neutron  induced 
changes  in  base  and  junction  voltages  at  high  injection  to  obtain 
theoretical  results  that  agree  with  the  current-voltage  characteristics 
of  actual  diodes,  including  neutron  effects.  With  a modest  history  of 
irradiation,  the  base  voltage  change  with  neutron  exposure  will  domi- 
nate, and: 

j • 7 ♦ K 0 (1) 

o 

with 

L = >/dT  (2) 


22.  J.M.  Swartz  and  M.O.  Thurston,  "Analysis  of  the  Effect  of  Fust- 
Neutron  Bombardment  on  the  Current-Voltage  Characteristic  of  a 
Conductivity-Modulated  p-i-n  Diode",  J.  Appl.  Phys.  , 37_,  No.  2, 
745-755,  (1966). 
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2kT  M . 1 W 

and  Vba»e  * V sinh  (3T}  tan  lsinh^) 


(3) 


whore 


t » initial  carrier  lifetime 
o 

t - carrier  lifetime  after  irradiation 
K ■ damage  constant 
f ■ neutron  fluence 
P ■ diffusion  constant 
L » ambipolar  diffusion  length 
k - Boltzman's  constant 
T ■ temperature 
e • electron  charge 

W » width  of  the  base  region  of  the  diode. 


Equations  1,  2 and  3 relate  the  observed  diode  property  (base  voltage], 
carrier  lifetime,  and  neutron  fluence.  These  relations,  and  the  life- 

1 8 

time  measurements  presented  by  Speers,  indicate  clearly  why  the  voltage 
observation  is  a measure  of  lifetime  degradation.  However,  the  explicit 
correctness  of  the  formulae  is  not  required  for  the  damage  measurements 
here.  With  the  allowance  of  interpolations  on  the  order  of  10\  of  the 
observed  changes,  damage  from  neutrons  from  a particular  source  is 
compared  with  an  equal  amount  of  damage  from  14  MeV  neutrons. 


The  SO-mil  diode  exposure-readout  procedure,  developed  for  measure- 
ments made  with  monoenergetic  neutrons,1*'  consisted  of  an  initial 

11  i 

exposure  to  14  MeV  neutrons  (about  10  n/cm*),  necessary  to  put  the 

diodes  in  a response  range  where  damage  measurements  are  the  result  of 
neutron  exposure  rather  than  material  and  device- fabrication  history. 
After  this  initial  irradiation,  the  diodes  were  annealed  for  several 
hours  at  a temperature  of  at  least  50°C  and  then  they  were  allowed  to 
stabilize  at  room  temperature  for  at  least  one  week.  At  a subsequent 
exposure,  the  diodes  were  irradiated  with  a fluence  of  neutrons  with  a 
selected  energy  spectrum  (fusion  or  fission  simulatorl.  Following  each 
such  exposure,  the  forward  voltage  was  measured  at  selected  times  for  a 
period  well  in  excess  of  100  hours  and  the  value  of  the  forward  voltage 
100  hours  after  mid-exposure  time  was  determined.  Typical  before  and 
after  voltages  were  1.2  and  1.4  volts.  The  individual  readings  were 
t 2 mV  and  thus  a 200  mV  damage- induced  change  could  be  determined  to 
about  t 2\  (this  does  not  include  the  errors  in  determining  neutron 
fluences).  The  diodes  were  exposed  in  sets  of  5 and  they  were  main- 
tained, except  for  the  period  of  transport,  at  30°C  after  exposure  to 
control  annealing.  A very  long-term  (8000  hours)  anneal  curve,  at  30#C. 
has  been  measured  for  these  devices.  This  curve  was  used  to  make  few- 
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percent  corrections  to  diode  voltages  when  they  were  returned  from 
exposure  too  late  to  read  at  the  100-hour  point.  The  nonlinear  response 
of  the  diodes  was  corrected  for  by  a polynomial  fit  to  experimental 
calibration  (14  MeV)  data,  and,  of  course,  this  nonlinearity  precludes 
quoting  a single  sensitivity  factor. 

The  neutron  fluence  measurement  techniques  include  use  of  acti- 
vation foils,  a proton-recoil  telescope,  and  an  associated-particle 
counter  at  the  BRl  Cockcroft-Walton  accelerator.  A weighted-average 
of  these  techniques  was  used  for  the  initial  evaluation  of  diode  response 

to  14  MeV  neutrons.  Also,  there  was  a subsequent  exposure  of  diodes  to 

23 

the  intense  14-MeV  neutron  source  at  Livermore.  For  equal  amounts  of 
damage,  the  fluences  for  the  Livermore  exposures  were  found  to  be  5 per- 
cent less  (0.947).  Because  it  is  a more  generally  available  source, 
all  fluences  were  corrected  to  agree  with  the  Livermore  value.  The 

fluence  measurements  at  the  APRF  reactor  were  obtained  using  techniques 

•>4 

given  by  McGarry  et  al.~  The  fluences  for  exposures  at  a cal ifornium- 
252  source  were  determined  by  exposure-time  and  position  with  a prior 
. . 24 

calibration  against  APRF.  In  all  these  measurements  advantage  was 

. 25 

taken  ot  very  recent  spectrum  evaluations  by  McGarry  et  al.  The 

importance  of  spectra  in  fluence  determinations  has  been  discussed  by 

26 

Dunn,  Kazi,  and  Saccenti.  The  fluence  measurements  for  exposures 
at  the  DORF  reactor  are  based  on  older  spectrum  measurements  and  may 
be  less  reliable  as  a result. 

Annealing  is  an  important  aspect  of  damage  measurement  and  three 
areas  of  concern  will  be  mentioned  with  regard  to  the  50-mil  measurements. 
First,  a stable  pre-irradiation  evaluation  and  a time  and  temperature 
controlled  post-irradiation  evaluation  are  required.  Second,  if  anneal- 
ing depended  on  neutron  energy,  then  the  energy  dependence  would  depend 

+ The  authors  would  like  to  thank  E.D.  McGarry  and  P.A.  Trimmer  of 
HDL  for  including  this  exposure  in  their  experiments. 

23.  R.  Booth,  H.H.  Barschall,  and  E.  Goldberg,  "Rotating  Target  for 
Intense  14-MeV  Neutron  Source",  IEEE  Trans.  Nucl.  Sci.,  NS-20, 

No.  3,  472-474,  June  (1973). 

tt  The  exposures  and  fluence  measurements  at  the  APRF  were  done  by  the 
reactor  staff  under  the  direction  of  Dr.  A.H.  Kazi. 

24.  E.D.  McGarry,  A.H.  Kazi,  G.S.  Davis,  and  D.M.  Gilliam,  "Absolute 
Neutron  Flux  Measurements  at  Fast  Pulse  Reactors  With  Calibration 
Against  Californium-252",  IEEE  Trans.  Nucl.  Sci.,  NS-23,  No.  b, 
2002-2006,  December  (1976). 

25.  E.D.  McGarry,  C.R.  Heimbach,  A.H.  Kazi,  and  G.W.  Morrison,  "Fast 
Pulse  Reactor  Neutron  Spectrum  Measurement  and  Calculations", 
preprint  (unpublished)  July  (1977). 

26.  T.A.  Dunn,  A.H.  Kazi,  and  J.  Saccenti,  "Fluence-to-Dose  Conversion 
Factors  for  APRF  Fast  Pulse  Reactor  Neutron  Spectra",  OSABRL  Report 
No.  1832,  September  (1975).  (AD  #B007605L) 


on  the  time.  We  have  proven  this  is  not  the  case.*  Third,  annealing 
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due  to  carrier  injection,  such  as  that  noted  by  Gregory  and  Sander, 

28  29 

by  Barnes,  and  by  Mai Ion  and  Harrity,  is  an  area  we  have  not  in- 

vestigated. It  can  be  said,  however,  that  a consistent  amount  of 
injection  was  involved  in  all  the  (50-rail)  measurements  presented. 

Very  great  care  is  necessary  in  characterizing  permanent  damage  in 
neutron  irradiated  devices.  The  stable  damage  observed  at  room  tempera- 
ture many  hours  after  irradiation  depends  in  a complex  way  on  the  time, 
temperature,  and  charge  injection  history.  As  indicated  in  the  para- 
graphs above,  for  the  50-mil  measurements  each  of  these  (time,  tempera- 
ture, and  injection  history)  was  kept  constant.  The  result  was  to 
require  a very  lengthy  and  careful  series  of  measurements.  To  provide 
an  alternate  procedure,  and  additional  measurements,  a different  test 
technique  was  used  with  a set  of  thirty-six  30-mil  diodes.  Rather  than 
maintain  a constant  temperature  in  an  oven,  the  30-mil  diodes  were  kept 
in  a controlled  lab  (24  to  28°C)  and  careful  measurements  of  the 
temperature  at  readout  were  used  to  make  small  (few-tenths-of-a-percent) 

30 

corrections.  A procedure  suggested  by  Kruger  et  al.  was  used  to  } 
eliminate  the  need  to  interpolate  to  a 100-hour-after- exposure  value. 
These  authors  found  (for  diodes  from  a different  manufacturer)  that 
diodes  annealed  for  2 minutes  at  100°C  gave  a constant  result  from  1 
to  300  hours  after  exposure.  The  procedure  used  was  to  expose  the 
diodes  to  100°C  (boiling  water,  the  diodes  were  kept  dry  by  plastic 
foil)  for  2.5  minutes  and  to  determine  the  forward  voltage  change  20 
days  after  exposure.  Tests  indicated  that  the  observed  forward  voltage 
change,  when  compared  with  the  value  20  days  after  irradiation,  was  2% 
high  at  10  days  and  2%  low  at  40  days.  Data  at  20  days  after  exposure 
was  available  for  all  runs  so  no  correction  was  required  for  time  of 
observation. 


27.  B.L.  Gregory  and  H.H.  Sander,  "Injection  Dependence  of  Transient 
Annealing  in  Neutron- Irradiated  Silicon  Devices",  IEEE  Trans. 

Nucl.  Sci.,  NS- 14 , No.  6,  116-126,  December  (1967). 

28.  C.E.  Barnes,  "Thermal  and  Injection  Annealing  of  Neutron- Irradiated 
P-type  Silicon  Between  76°K  and  300°K",  IEEE  Trans.  Nucl.  Sci., 
NS-16,  No  6,  28-32,  December  (1969). 

29.  C.E.  Malion  and  J.W.  Harrity,  "Short-Term  Annealing  in  Transistors 
Irradiated  in  the  Biased-off  Mode",  IEEE  Trans.  Nucl.  Sci.,  NS-18, 
No.  6,  45-49,  December  (1971). 

30.  H.  Kruger,  G.  Tumbragel,  R.  Metzner,  and  H.  Koch,  "Fast-Neutron 
Dosimetry  With  Silicon  Diodes",  Neutron  Monitoring  for  Radiation 
Protection  Purposes,  Vol.  II,  IAEA-SM-167/53,  401-409,  Inter- 
national  Atomic  Energy  Agency,  Vienna  (1973). 
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The  injection  dependence  of  the  time,  temperature,  Injection 
history  influence  on  annealing  is  the  least  understood  aspect  of  this 

complex  interaction.  In  addition  to  the  studies  previously  noted,*  ’ * 

an  increase  in  forward  voltage  or  "reverse  anneal"  has  been  observed.  1 
An  investigation  of  this  effect  showed  that  the  forward  voltage  after 
irradiation  land  thus  the  measured  voltage  change)  will  increase 
significantly  with  continued  current  injection  for  protracted  periods, 
for  example,  several  diodes  with  very-permanent  neutron  damage  UJ 
years  after  irradiation)  were  found  to  increase  from  1.19  to  1 . -S  V 
when  injected  for  several  hours.  This  is  a 51  increase  in  the  forward 
voltage  and  a 151  increase  in  tho  neutron- induced  voltage  change,  lhis 
effect  was  found  to  saturate  after  150  to  300  minutes  of  injection,  and 
it  was  found  to  increase  with  an  increase  in  the  amount  of  neutron 
damage.  No  effect  was  observed  for  unirradiated  diodes  and  the  effect 
was  reduced  when  some  neutron  damage  was  annealed  at  high  temperature. 

A few  months  at  room  temperature,  or  a few  minutes  at  100"C,  was  found 
to  bleach  this  effect.  The  30-mil  diodes  were  Injected  (standard 
forward  current  of  0.1A)  for  a period  in  excess  of  300  minutes  in  sets 
of  10.  They  were  then  transferred  to  a different  circuit  conf igurat ion 
and  read  at  5,  30,  and  50  seconds  after  the  start  of  current  Injection. 
The  30s  values  were  used  to  determine  damage  and  the  other  values  were 
used  to  insure  that  injection  effects  were  stabilized. 

To  summarize  the  diode  exposure- readout  procedures,  the  50  mil 
diode  measurements  were  made  with  the  same  techniques  used  for  studies 
with  monoenerget ic  neutrons.  These  techniques  are  described  in  the 
first  paragraphs  of  this  chapter,  and  they  have  been  presented  in 

previous  reports.**'  The  same  careful  procedures  were  followed  with 
the  30-mil  diodes.  However,  the  30-mil  diodes  were  read  after  they 
had  been  annealed  J.5  minutes  at  100*1),  after  they  had  been  current 
injected  for  at  least  300  minutes,  and  after  JO  days  at  room  tempera 
t lire  subsequent  to  exposure,  Exact  readout  temperature  and  inieetton 
dependent  stability  were  noted. 


in.  cai.cumtionai.  m>i:i;i>imi;s 

In  our  first  efforts  to  determine  the  energy  dependence  of  neutron 

l IJ 

induced  displacement  damage  in  silicon  *'  a limited  (unevaluated)  set 
of  cross  section  data  for  specific  neutron  reactions  was  used  with  an 


31.  T.  (Trober , Private  Communication. 

^J.  W.R.  Van  Antwerp  and  d.E.  Youngblood,  "displacement  Damage  in 
Silicon  Irradiated  With  n-  to  10-MeV  Neutrons",  Bull.  Am.  Plus. 
Soc.,  1^*,  p.  o'* 8.  dune  ( 1 972) . 
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approximation  for  the  LSS  displacement  fraction  taken  from  Bertin 
xx  9 

et  al.  This  approximation,  an  empirical  fit  to  the  LSS  nuclear 
energy  loss,  is  given  by: 


dc 
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nucl 
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0. 67* 2. 07c *0. 03c 


(4) 


where  c and  p are  the  dimensionless  parameters  defined  by  LSS  in 
terms  of  particle  energy,  E,  and  displacement,  X,  as 


with 
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Z^e^A^A,)  (A^A2r 

Z^.Aj  - atomic  number  and  mass  of  projectile 
Z2,A , - atomic  number  and  mass  of  stopping  medium 
N - number  of  stopping  atoms  per  unit  volume 

a ■ 


0.8853  a / (Z.  2/3«-Z,2/3)lj 

Oi  a 
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a^  - the  Bohr  radius,  0.529  X 10  cm 
e ■ the  electron  charge. 


(5) 


The  LSS  electronic  energy  loss  was  taken  as  proportional  to  the  square 
root  of  the  energy  parameter,  c,  and  the  formulation  of  an  integral  for 
the  nuclear-stopping  fraction  followed  in  a straightforward  way  (details 
can  be  found  in  reference  2). 

In  integrating  the  formulae  for  energy,  nuclear  stopping  is  being 
treated  as  a continuous  process  and  this  is  an  approximation  as  the 
energy  lost  in  a single  collision  can  be  a substantial  fraction  of 
the  total  recoil  energy.  Also,  LSS  as  used  in  the  preliminary  cal- 
culations provides  the  fraction  of  onergy  of  the  primary  recoil  that 
goes  into  nuclear  scattering  and  no  correction  was  made  in  the  pre- 
liminary calculations  for  multiple  collisions.  That  is,  no  correction 
was  made  for  energy  lost  to  ionization  in  secondary,  tertiary,  et  seq. , 
collisions.  The  preliminary  program  was  designed  with  particular 
emphasis  on  determining  the  effects  of  various  parts  of  the  input 
neutron  cross  section  information  and  it  served  this  purpose  well. 
Despite  its  shortcomings,  it  also  served  as  a base  for  the  more 
elaborate  subsequent  effort. 


3T!  M.C.  Bertin,  N.  Benczer-Koller,  G.G.  Seaman,  and  J.R.  MacDonald, 

Cg 

"Electromagnetic  Transition  Rates  in  ' Ni",  Phys.  Rev.,  183, 
9M-977,  (1909). 
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The  final  program  for  damage  calculations  was  formulated  with  a 
desire  to  improve  the  evaluation  of  the  displacement  fraction  and  with 
a continued  interest  in  the  sensitivity  of  damage  calculations  to  cross 
section  inputs.  The  ability  to  read  an  evaluated  cross  section  set 
from  magnetic  tape  was  considered  necessary,  and  there  was  a need  for 
very  good  neutron  energy  resolution  for  use  in  analyzing  accelerator 
neutron  damage  experiments. 

The  partial  cross  section,  solid  angle,  and  primary  recoil  energy 
calculations  were  done  in  the  same  way  as  for  the  preliminary  calculations, 
but  with  the  neutron  cross  sections  and  angular  distributions  taken  from 
a tape.  The  preliminary  program  calculated  the  final  displacement  energy 
for  each  neutron  reaction  at  each  of  60  angles  and  summed  the  results. 

The  present  program  made  similar  calculations,  then  doubled  the  number 
of  angles  and  repeated  the  calculation  with  a test  (0.5%)  for  conver- 
gence, with  further  doubling  of  the  number  of  angles  when  required. 

Also,  there  was  a substantive  change  in  the  calculation  of  the  frac- 
tion of  the  primary  recoil  energy  available  for  displacements.  A 
correction  was  made  for  the  energy  lost  to  ionization  in  the  collision 
sequence  and,  for  this,  a numerical-integration  technique  was  used  in 
an  iterative  process.  Many  collisions  were  followed  for  the  higher- 
energy  recoils  and  several  tens  of  values  were  determined,  geometri- 
cally spread  across  each  decade  of  recoil  energy.  A few  hundred  values 
were  stored,  then,  for  table  lookup  and  interpolation. 

Calculations  of  the  silicon  recoils  in  silicon  were  supplemented 
with  calculations  of  the  aluminum  and  magnesium  recoils  in  silicon  to 
account  for  the  heavy  recoil  products  of  the  (n,ci)  and  (n,p)  reactions. 

The  proton  and  alpha  recoils  were  not  included  in  the  final  calculations. 
There  are  two  reasons,  other  than  convenience  and  simplicity,  for  omitting 
these  light-particle  recoils.  First,  the  contributions  from  these  light 
particles  were  evaluated  in  a few  cases  and  found  to  be  no  more  than  1 or 

2%  of  the  total  damage  (similar  results  were  obtained  by  Smith  et  al.1)- 
Second,  these  particles  have  ranges  significant  when  compared  to  typical 
device  dimensions,  and  a realistic  calculation  of  the  p and  a damage 
wouid  have  to  account  for  losses  and  gains  of  p's  and  a's  with  the  surround 
ing  materials.  The  cross  section  tape  used  was  the  DNA  file  designated 
as  DNA  4151  MOD  3.  Except  for  editorial  changes,  this  cross  section 
information  is  the  same  as  that  in  ENDF/B,  Version  IV,  MAT  1194.  The 
calculations  include  elastic  scattering,  inelastic  scattering  (22 
excited  states  continuum),  (n,p)  reaction  (15  excited  states+continuum) , 
(n,a)  reaction  (12  excited  states+continuum) , (n,2n),  (n,np),  and 
(n,n«).  In  the  last  three  reactions,  the  heavy  reaction  product  was 
assumed  to  be  at  rest  in  the  center  of  mass  frame,  and,  in  all  cases, 
the  reaction  was  assumed  to  be  isotropic  in  the  center  of  mass  system 
if  no  other  information  was  available. 


to 


Although  the  above  'present-work'  formulation  provided  excellent 

34 

results,  a simple  approximation  taken  from  Doran'  will  produce  almost 
identical  results,  and  it  seems  most  important  to  point  out  that  this 
part  of  the  computation  does  not  contribute  to  differences  between  the 
present  results  and  those  of  others.  The  formula  given  by  Doran,  when 
reduced  to  a single  constituent  (silicon  recoil  in  silicon)  is: 


f l*kg(Er) 

(6) 

with 

k - 0.1334Z2/V1/2 

(7) 

and 

g(E  ) - E +0.40244E  3/4  + 3.4008E  1//f> 
r r r r 

(8) 

where  f ■ fraction  of  recoil  energy  going  into  displacements, 
corrected  for  multiple  collisions. 

Z,A  * atomic  number  and  mass  of  silicon  (14,28). 

Er  = recoil  (or  PKA,  primary  knock-on  atom)  energy. 

Table  I provides  a comparison  of  the  displacement-fraction  estimates 
obtained  from  Equation  (4),  from  present-work,  and  from  Equation  (6). 


TABLE  I.  Comparison  of  Estimates  of  Displacement  Energy  F.j , at  Several 
Recoil  Energies,  Ef  (all  values  in  keV) . 


E 

r 

(keV) 

Ed 

(Equation  4) 

Ed 

(Present -Work) 

Ed 

(From  Equation  6) 

1 

0.91 

0.79 

0.78 

3 

2.7 

2.3 

2.24 

10 

8.8 

7.0 

6.9 

30 

25. 

19. 

18.4 

100 

68. 

50. 

49. 

300 

141. 

101. 

99. 

1000 

240. 

167. 

165. 

3T!  b.C.  Doran,  "Neutron  Displacement  Cross  Sections  for  Stainless 
Steel  and  Tantalum  Based  on  a Lindhard  Model",  Nucl.  Sci.  Eng., 
49,  No.  2 130-144,  October  (1972). 
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The  table  shows  both  the  importance  of  the  multiple-collision  correc- 
tion and  the  close  correspondence  of  the  present  calculations  to  those 

34 

(including  any  future  ones)  baaed  on  the  Doran  approximation.  This 
approximation,  Equation  (6),  is  in  a convenient  form  and  equally  con- 
venient expressions  could  be  obtained  for  aluminum  recoils  in  silicon 

and  magnesium  recoils  in  silicon.  The  comparison  between  present-work 

•> 

and  Doran  is  somewhat  better  than  the  previous  used  comparison'  taken 
35 

from  Mayer  et  al. 

The  results  of  the  calculations  on  energy  dependence,  in  complete 
tabular  form,  are  required  in  order  to  calculate  the  displacement 
damage  to  be  expected  from  an  accelerator  produced  neutron  experiment, 
from  a simulator  spectrum,  or  from  a threat  spectrum,  for  these  uses 

2 

the  tables  of  calculated  results,  presented  in  a previous  report,'  are 
included  in  the  appendix.  The  calculations  were  done  for  each  5 keV 
interval  from  0.1  to  20  MeV  and  these  results,  for  0.1  to  16  MeV,  are 
given  in  the  first  four  tables.  There  are  also  tables  providing  25 
keV  interval  averages  for  0.1  to  16  MeV  and  100  keV  interval  average 
values  for  the  energy  range  0.1  to  20  MeV.  The  appendix  also  contains 
tables  of  the  calculated  response  for  various  group  structures  used  in 
subsequent  sections  of  this  report. 


IV.  COMPARISON  OF  CALCULATIONS  AND  MEASUREMENTS 

Prior  to  presenting  the  results  for  more  complex  neutron  spectra, 
some  of  the  calculations  and  measurements  done  for  accelerator-produced 
neutrons  will  be  compared.  This  will  provide  the  reader  with  a better 
feel  for  the  level  of  agreement  that  might  be  expected  between  calcu- 
lated and  observed  values.  A graphic  presentation  of  the  totality  of 
calculated  results  is  given  in  Figures  1 and  2.  The  fractional  contri- 
butions from  the  various  neutron  reactions  are  shown  as  well  as  the 
total  calculated  damage.  Figure  1 covers  the  neutron  energy  range 
0.1  to  10  MeV  while  Figure  2 covers  10  to  20  MeV.  To  calculate  the 
expected  damage  the  (S-keV  interval)  tabulated  damage  values  are 
averaged  over  the  appropriate  energy  interval  (for  the  energy  range  of 
an  accelerator  run)  or  over  each  of  the  appropriate  energy  intervals 
(for  a given  group  structure).  The  resulting  damage  factor (s)  is 
multiplied  by  the  neutron  fluence  (or  fraction  of  a neutron  in  each 
group) . 


35 . J.W.  Mayer,  L.  Eriksson,  and  J.A.  Davies,  Ion  Implantation  in 
Semiconductors , p.  68,  Academic  Press,  NY,  (1970). 
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Figure 


The  results  of  damage  calculations  are,  typically,  0.1  to  0.2 
MeV*b  (100  to  200  MeV*m  b) , with  the  specific  value  of  0.187  MeV*b 
for  14.2  MeV  neutrons.  The  techniques  used  to  process  experimental 

data,*’"*  corrected  for  non-linear  response  and  normalized  to  a neutron 

10  2 

fluence  of  2.357  X 10  m/cm  , produced  measured  damage  results  that 
ranged  from  1.4  to  5.7.  The  measured  values  divided  by  calculated 
values  for  the  same  neutron  energy  (e.g.,  5/0.2)  produced  numbers 
of  the  magnitude  25,  and  the  constancy  of  this  ratio  is  an  indication 
of  the  adequacy  of  the  model  used  to  calculate  the  damage.  A comparison 
of  calculated  and  measured  results  for  neutron  energies  near  1 MeV  is 
illustrative  of  the  facts  just  described  and  these  are  shown  in  Table  II. 


TABLE  II.  Comparison  of  Calculated  and  Measured  Damage  Near  1 MeV. 


Neutron  Energy 
(MeV) 

3 ) 

Energy  Spread  J 
(keV) 

Calc.  Dam. 
(MeV*b) 

Ratio*3) 

(Exp/calc) 

Error 

(%) 

0.696 

60 

0.0582 

24.4 

13 

0.957 

49 

0.114 

23.7 

13 

1.157 

44 

0.0502 

30.3 

13 

1.630 

38 

0.193 

26.4 

13 

2.370 

37 

0.108 

24.7 

13 

3.990 

35 

0.131 

20.8 

12 

4.990 

30 

0.152 

18.9 

12 

14.2 

- 

0.187 

22.7 

7 

/ 

j1) 

JThe  energy  spread  includes  line-shape  effects.  The  stated  energy 
widths  define  an  interval  that  contains  90%  of  the  total  neutrons 
incident  on  the  diode  samples. 


The  ratios  given  include  experimental  values  (see  text  for  units) 
obtained  using  fluences  measured  with  a long  counter  except  for 
14.2  MeV.  The  14.2  MeV  value  was  obtained  using  a fluence  determined 
by  proton-recoil  telescope. 

These  results  are  shown  graphically  in  Figure  3,  and  the  results  for 
4.0  to  7.5  MeV  for  experiments  done  at  the  BRL  tandem  Van  de  Graaff 
are  shown  in  Figure  4.  For  both  Figures  3 and  4 the  calculated  (left) 
and  experimental  (right)  scales  are  for  a ratio  of  26.1/1.  This  is 
the  average  value  for  the  25  results  shown. 
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Figure  4.  Calculated  and  Measured  Damage,  4 to  7.5  MeV 


It  has  been  common  practice  to  compare  all  experimental  or 
calculated  results  on  energy  dependence  of  neutron  damage  in  silicon 

to  a Holmes  curve1*  or  an  expression  such  as  that  suggested  by 

Messenger."5**  This  relation  is  of  the  form 

PIE)  » A(l-e'B/E)  19) 

where  0(E)  = displacement  damage 
E = neutron  energy 
A,B  = constants. 

Various  values  of  the  constants  A and  B have  been  chosen  to  fit  past 
theoretical  or  experimental  results.  Using  the  values  (from  Messenger) 
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quoted  by  Green  and  Thatcher,  A»l.l  and  B=2.2,  a Messenger  curve  is 
shown  in  Figure  5 overlaid  on  the  present  calculated  results  averaged 
over  100  keV  intervals.  The  Messenger  curve  has  been  normalized  to  a 
value  of  187  MeV • mb  at  14.2  MeV. 


V.  RESULTS 

Experimental  evaluations  of  the  fusion/fission  damage  ratios  for 
neut ron- induced  displacement  damage  in  silicon  have  been  done  with  the 

'S ' . 

Al’Rl  reactor  (both  glory  hole  and  leakage  spectral,  with  a "‘"ft 
spontaneous-fission  source,  and  with  the  PORE  reactor.  In  each  case 
a set  of  S,  1.27  mm  base-width,  diodes  and  a set  of  S,  0.  7b  nun  base- 
width,  diodes  were  each  exposed  to  a suitable  fluence.  The  techniques 
by  which  the  fluences  were  determined  and  the  two  methods  used  for 
evaluating  the  diode  responses  have  been  described  in  earlier  sections 
of  this  report.  The  results  of  these  icasurements  are  given  in  Table  111. 
Although  neutrons  in  the  energy  range  10  to  100  keV  do  not  create  sig 
nificant  amounts  of  displacement  damage  (at  least  not  for  threat  and 
simulator  spectra  considered  here),  the  fission  simulator  fluence 
measurements  are  for  the  number  of  neutrons  with  an  energy  greater  than 
10  keV  and  thus  10  keV  (or  9 keV)  calculated  values  should  be  used  for 
comparison.  A happen-stance  of  exposure,  near  the  maximum  of  a 
suggested  range,  made  it  necessary  to  extrapolate  the  nonlinear  diode- 
's ' 

response  curve  (by  a tew  percent)  to  the  value  tor  Ct"  with  0."o  mm 
base  -wide  diodes  and  this  result  should  be  given  slightly  less  weight. 


3b.  G.C.  Messenger,  "Displacement  Damage  in  Silicon  and  Germanium 

Transistors",  IEEE  Trans.  Nucl.  Sci.,  12,  No.  2,  5374  April  (I9b5). 
M.L.  Green  and  R.K.  Thatcher,  "Preparation  of  a Standard  Technique 
for  Determination  of  Neutron  Equivalence  for  Bulk  Damage  in 
Silicon",  IEEE  Trans.  Nucl.  Sci.,  NS-19,  No.  n,  200-208,  December, 
(1972). 
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DAMAGE  (CALCULATED,  MeV-mb) 


Figure  5.  Calculated  Damage  Averaged  Over  100  keV  Intervals 


FABLE  111.  Measured  Fusion/Fission  Damage  Ratios. 


Source 

Diode  Base 

Width  (mm) 

Damage  Rat io 
(14  MeV/souree) 

APRF  (Glory  Hole) 

1.27 

1.75 

II  It  M 

0.76 

1.84 

" (leakage) 

1.27 

1 . nt> 

H II 

0 . 7t> 

1.70 

:5Jcf 

1.27 

1.52 

II 

0.7b 

1.80 

PORF  (2cm  11,0) 

1.27 

2.45 

II 

0.7b 

2.58 

Calculated  damage  for  simulator  spectra  and  the  ratios  of  these  to 
the  value  at  14  MeV  have  been  obtained.  The  spectra  used  for  APRF 

•>5  •> 

(.(.'.lory  Hole),  APRF  (Leakage),  and  "Cf  were  94-group  spectra  supplied 
by  F.P.  McGarrv.  The  damage  factors  for  these  energy  groups  are  given 
in  the  appendix,  Table  A-X.  The  spectra  used  for  the  PORF  reactor 
calculations  were  taken  from  "TREE  Preferred  Procedures",  DNA  2028H. 
page  5-49,  19'2.  The  damage  factors  for  these  groups  are  given  in  the 
appendix,  Fable  A- XII.  Damage  ratios  for  source-neutron  energies,  T , 

greater  than  100  keV  and  greater  than  9 keV  are  given  in  Table  IV. 

Calculated  damage  for  several  threat  spectra  have  been  obtained. 

These  spectra"'5'  were  supplied  in  a .V*- group  energy  structure  with  the 
energy  intervals  given  in  the  appendix,  Table  A-Xl.  The  damage  factors 
for  these  groups  are  also  given  in  TABLE  A-XI.  These  spectra  are,  in 
each  case,  a neutron  spectrum  transported  over  an  air/ground  interface 
due  to  a fission  weapon  (Spectrum  A),  a tactical  nuclear  weapon  (Spec- 
trum R) , and  a modified  14  MeV  source  (Spectrum  C) . These  results, 
including  the  ratios  of  1-14  MeV  neutron  to  1-source  neutron  of  energy, 

F , greater  than  100  keV  and  greater  than  9 keV,  are  given  in  Table  \. 
aVso,  the  numbers  of  neutrons  in  energy- groups  down  to  sub-thermal 
were  given  so  damage  ratios  for  neutrons  of  all  energies  were  determined 
(column  4,  Table  V). 


38.  A.E.  Rainis  (Private  Communication). 


TABLE  IV.  Calculated  Ratios  of  Damage  for  Simulator  Spectra. 


Source 


Damage  Ratio 
(1-14  MeV  neutron  to 
1 source  neutron  with 
E > 100  keV) 


Damage  Ratio 

(1-14  MeV  neutron  to 

1 source  neutron  with 

E > 9 keV) 
n 


APRF  (Glory  Hole) 

2.13 

2.19 

APRF  (Leakage) 

2.07 

2.12 

2S2Cf 

1.80 

1.82 

DORF  (Center  of  Core) 

2.13 

2. 54 

DORF  (2  cm  in  H20) 

2.10 

2.54 

14  MeV 

1.00 

1.00 

TABLE  V.  Calculated 

Ratios  of  Damage  for  Typical  Threat 

Spectra. 

Damage  Ratio 

Damage  Ratio 

Damage  Ratio 

(14  MeV/ source) 

(14  MeV/source) 

(14  MeV/source) 

Source 

E > 100  keV 
n 

E > 9 keV 
n 

E > 0 keV 
n 

Threat  Spectrum  A 

2.58 

7.06 

10.4 

Threat  Spectrum  B 

2.74 

6.47 

9.40 

Threat  Spectrum  C 

1.98 

3.90 

5.32 

14  MeV 

1.00 

1.00 

1.00 

Spectrum  A = fission  weapon,  transported  over  air/ground  interface 
Spectrum  B = tactical  nuclear  weapon,  transported  over  air/ground 
interface 

Spectrum  C = modified  14  MeV  source,  transported  over  air/ground 
interface. 


VI.  DISCUSSION 

The  nature  of  modeling  permanent  neutron  damage  in  materials  can 
readily  be  divided  into  four  phases: 

a.  Model  the  nuclear  reaction  cross  sections,  combine 
these  with  measured  cross  sections,  and  form  a self- 
consistent  evaluated  neutron  cross  section  set. 
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b.  Apply  the  known  kinematics  (and  model  where 
necessary)  for  each  nuclear  reaction  at  each 
angle  and  obtain  the  energy  distribution  of 
primary  recoils. 

e.  Model  the  fractioning  of  primary-recoil  energy 
into  electronic  stopping  (not  producing  per- 
manent damage)  and  energy  available  for 
displacements. 

d.  Model  the  conversion  of  displacement  energy 

into  the  number  and  geometry  of  displaced  atoms. 


In  the  frame  of  these  four  phases,  the  calculations  of  damage  presented 
in  this  report  are  based  on  a particularly  thorough  and  detailed  con- 
sideration of  phases  a and  b.  These  phases  are  discussed  in  detail  in 
the  section  on  calculat ional  procedures.  Phase  c has  been  done  using 

9 

the  fractioning  ot  LSS,  and,  in  addition  to  the  confirmation  by 

10  39  40 

Sattler  and  Vook,  the  work  of  Chaseman,  et  al.,“'  has  shown  that 

the  ionization  losses  for  germanium  atoms  moving  in  germanium  agree 

with  LSS  (electronic  stopping  proportional  to  velocity)  for  the  energy 


range  10  to  100  keV 


Jones  and  Kraner 
4 


41 

have  extended  tins  with  measure 
with  a final  check  at  254  eV . liven 


ments  in  the  range  1 to  2 keV  and 

4 2 

at  254  eV  Tones  and  Kraner  “ find  an  inelastic  loss  of  ca.  40  eV  in 
good  agreement  with  LSS.  With  regard  to  phase  d,  it  has  been  assumed 
that  the  number  of  defects  initially  formed  is  proportional  to  the 
energy  available  for  displacements,  that  the  number  of  permanent  de- 
fects is  proportional  to  the  number  initially  formed  and  independent 
of  their  geometry,  and  that  the  mechanical  or  electrical  property  of 
interest  does  not  depend  on  the  geometry  (clustering)  of  defects. 


T9,  G.  Ghasman,  K.W.  Tones,  and  R.A.  Ristinen,  "Measurement  of  the 
Lnergy  Loss  of  Germanium  Atoms  to  Electrons  in  Germanium  at 
Lnergies  Below  100  keV,"  Phys.  Rev.  Letters,  1^5,  No.  t>,  245-249, 
August,  (19t>5). 

40.  G.  Ghasman,  K.W.  Tones,  and  H.W.  Kraner,  "Band-Gap  P.ffects  in 

7 ->* 

the  Stopping  of  Ge' “ Atoms  in  Germanium",  Phys.  Rev.  Letters, 

2J,  No.  20,  1430-1433,  November  (19o8).  ^ 

41.  K.W.  Tones  and  H.W.  Kraner,  "Stopping  of  1-  to  1.8-keV  ' ' Ge  Atoms 
in  Germanium",  Phys.  Rev.  G,  4,  125-129,  Tanuary,  (1971). 

42.  K.W.  Tones  and  H.W.  Kraner,  "Lnergy  l.ost  to  Ionization  by  254  eV 

' Ge  Atoms  Stopping  in  Ge",  Phys.  Rev.  A,  11,  No.  4,  1347-1353, 
Apri 1 , (1975) . 


The  measured  values  of  damage  ratios  are  primarily  limited  by 
fluence  determinations.  Although  neutron- induced  voltage  changes  can 
be  measured  to  1 or  2 percent,  an  error  of  10\  is  possible  in  the 
fluence  evaluations  at  14  MeV.  For  pulsed  reactor  neutron  fluence 
measurements  the  neutron  spectrum  is  involved.  The  spectra  are  used, 
for  example,  to  determine  the  fluence  ratio,  U(En  > 10  keV)/0(E  > 3 
MeV)  and  this  ratio  is  used  to  convert  activation  foil  measurements 
to  fluences. 

The  experimental  results  in  Table  III  show  slight  but  consis- 
tently higher  damage  ratios  from  the  0.76  mm  diodes.  If  the  value 
252 

for  Cf  is  discounted,  the  difference  is  about  4\  and  cannot  be 
considered  significant.  This  would  be  worth  a careful  check  if  the 
exposures  are  repeated,  however.  When  the  measurements  in  Table  III 
are  compared  with  the  calculated  (E  > 9 keV)  results  in  Table  IV, 
spectrum  and  fluence  measurements  aDe  again  brought  into  question. 

252 

The  three  results  on  the  APRF  and  ‘ Cf  simulators  are  based  on  the 
recent,  carefully  determined  94-group  spectra  and  they  are  self-con- 
sistent to  about  7\,  hut  the  measured  ratios  are  15l  less  than  the 
calculated  ratios.  The  more  limited  and  dated  spectral  information  on 
the  DORF,  meanwhile,  has  produced  an  experimental-to-calculated  ratio 
of  1.00.  It  has  been  noted  that  (especially  with  non-monoenerget ic 
sources)  the  spectrum  used  affects  both  the  calculated  damage  and  the 
fluence  measurement . 


VII.  CONCLUSIONS 

Despite  the  limitations  cited  in  the  discussion  above,  excellent 
results  have  been  obtained  on  the  energy  dependence  of  neutron  damage 
in  silicon  and  a modeling  technique  of  general  applicability  to  damage 
calculations  has  been  presented.  Neutron- induced  displacement  damage  in 
silicon  has  been  calculated  and  measured  for  a representative  selection 
of  simulator  and  threat  neutron  spectra.  The  calculated  damage  is  based 
on  the  Lindhard  nuclear  stopping  fraction  corrected  for  multiple  col- 
lisions and  it  thus  provides  the  energy  available  for  displacement 
damage.  Tables  of  calculated  damage  in  5,  25  and  100  keV  intervals  are 
presented  as  well  as  tahles  of  calculated  damage  for  several  frequently 
used  group  structures.  An  experimental  check  of  the  damage  has  been 
obtained  using  wide-base  silicon  diodes  at  a high  injection  level. 

The  changes  in  forward  voltage  at  a fixed  current,  reflecting  the 
degradation  of  injected  carrier  lifetime,  were  observed. 

The  uncertainties  in  the  measurements  we  have  presented  of  fusion/ 
fission  damage  ratios  for  neutron  induced  displacement  damage  in 
silicon  are  generally  less  than  20V  The  comparison  of  measurements 
and  calculations  is  also  less  than  20\  and  this  is  sufficient  for  most 
military  applications.  The  problems  with  a 1 MeV  hardness-assurance 


criteria,  delineated  most  concisely  in  an  earlier  report,  remain. 
However,  if  adequate  spectral  information  is  folded  into  the  damage 
table  presented,  any  equivalence  definition  can  be  accommodated. 
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TABLES 


The  basic  damage  calculations  were  done  in  5-keV  intervals  from 
0.1  to  20  MeV  and  the  results,  stored  on  disc,  remain  available.  To 
provide  the  information  necessary  for  others  to  obtain  a set  of  damage 
factors  for  any  neutron  energy  group  structure  of  interest,  the  complete 
S-keV  interval  results  from  0.1  to  16  MeV  are  presented  in  Tables  A-I, 
A-II,  A- III  and  A-IV.  Tables  A-V  and  A-VI  provide  25-keV  interval 
averages  for  0.1  to  16  MeV,  and  Table  A-VII  has  100-keV  average  values 
for  0.1  to  20  MeV.  The  results  for  a 22-group  structure  (Table  VIII), 
a 94-group  structure  (Table  IX),  a 36-group  structure  (Table  X)  and  a 
24-group  structure  (Table  XI)  are  also  presented. 


TABLE  A-l.  Relative  Neutron  Damage  (MeV-b)  Averaged  Over  S keV  Intervals 
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TABLE  A-II.  Relative  Neutron  Daaage  (MeV.  b)  Averaged  Over  5 keV  Intervals  - 4 to  8 MeV 
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TABLE  A- VI II. 

Calculated  Damage  Results  in 

22-Group  Structure. 

Group 

AE 

n 

Damage 

Number 

(MeV) 

(MeV* mb) 

1-9 

E <•  1 
n 

-- 

10 

.11  - .55 

53 

11 

.55  - 1.11 

84 

12 

1.11  - 1.83 

97 

15 

1.83  - 2.35 

125 

14 

2.35  - 2.4b 

103 

15 

2.46  - 3.01 

127 

10 

3.01  - 4.07 

120 

17 

4.07  - 4.97 

151 

18 

4.97  - b. 3b 

153 

19 

b . 3b  - 8.19 

158 

20 

8.19  - 10. 0 

16b 

21 

10.0  - 12.2 

17b 

12.2 


15.0 


18? 


TABLE  A-IX.  Calculated  Damage  Results  in  94-Group  Structure. 


Group 

Number 

*En 

(MeV) 

Damage 

(MeV • mb) 

1-25 

V1 

-- 

26 

.100-. Ill 

5 

27 

.111-. 136 

4 

28 

.136-. 183 

18 

29 

. 183- . 200 

104 

30 

.200-. 224 

91 

31 

.224-. 250 

67 

32 

.250-. 302 

55 

33 

.302-. 400 

51 

34 

.400-. 408 

52 

35 

.408-. 498 

54 

36 

.498-. 600 

80 

37 

.600-. 672 

54 

38 

.672-. 800 

66 

39 

.800-. 821 

165 

40 

.821-1.00 

101 

41 

1.00-1.11 

77 

42 

1.11-1.20 

55 

43 

1.20-1.35 

81 

44 

1.35-1.40 

82 

45 

1.40-1.50 

103 

46 

1.50-1.60 

106 

47 

1.6 0-1.80 

124 

48 

1.80-2.00 

151 

49 

2.00-2.20 

107 

50 

2.2 0-2.23 

108 

51 

2.23-2.30 

108 

52 

2.30-2.40 

108 

53 

2.40-2.60 

131 

54 

2.60-2.73 

124 

55 

2.73-2.80 

111 

56 

2.80-3.00 

122 

57 

3.0 0-3.20 

126 

58 

3.20-3.40 

118 
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TABLE  A- IX 


(Continued) 


Group 

Number 

AE 

n 

(MeV) 

Damage 

(MeV -mb) 

59 

3.40-3.60 

120 

60 

3.60-3.68 

78 

61 

3.68-3.70 

84 

62 

3.70-3.80 

119 

63 

3.80-4.00 

125 

64 

4.00-4.20 

125 

65 

4.20-440 

152 

66 

4.40-4.60 

145 

67 

4.60-4.80 

172 

68 

4.80-5.00 

158 

69 

5.00-5.20 

162 

70 

5.20-5.40 

135 

71 

5.40-5.60 

136 

72 

5.60-5.80 

167 

73 

5.80-6.00 

155 

74 

6.00-6.07 

161 

75 

6.07-6.50 

153 

76 

6.50-7.00 

147 

77 

7.00-7.50 

160 

78 

7.50-8.00 

172 

79 

8.00-8.50 

161 

80 

8.50-9.00 

165 

81 

9.00-9.50 

166 

82 

9.50-10.0 

168 

83 

10.0-10.5 

172 

84 

10.5-11.0 

172 

85 

11.0-11.5 

175 

86 

11. 5-12.0 

183 

87 

12.0-13.0 

186 

88 

13.0-14.0 

188 

89 

14.0-14.9 

187 

90 

14.9-15.0 

187 

91 

15.0-16.0 

187 

92 

16.0-17.0 

186 

93 

17.0-18.0 

183 

94 

18.0-20.0 

184 
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TABLE  A-X.  Calculated  Damage  Results  in  37-Group  Structure. 


droup 

Number 

AEn 

(MeV) 

Damage 

(MeV*mb) 

1-14 

E <•! 
n 

-- 

15 

.111-. 158 

4 

16 

.158-. 550 

59 

17 

.SS0-1.00 

85 

18 

1.11-1.83 

96 

19 

1.83-2.31 

127 

20 

2.31-2.39 

109 

21 

2.39-3.01 

124 

22 

3.01-4.07 

120 

23 

4.07-4.72 

144 

24 

4.72-4.97 

170 

25 

4.97-6.38 

153 

26 

6.38-7.41 

152 

27 

7.41-8.19 

168 

28 

8.19-9.05 

166 

29 

9.05-10.0 

166 

30 

10.0-11.1 

172 

31 

11.1-12.2 

180 

32 

12.2-12.8 

186 

33 

12.8-13.8 

188 

34 

13.8-14.2 

188 

35 

14.2-14.9 

187 

36 

14.9-16.9 

187 

37 

16.9-19.6 

184 
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TABLE  A'XI . Calculated  Damage  Results  in  24-Group  Structure. 


1 


Group 

Number 

AE 

n 

(MeV) 

Damage 

(MeV* mb) 

1 

10.0  - 7.79 

166 

") 

7.79  - 6.07 

156 

3 

6.07  - 4.72 

152 

4 

4.72  - 3.68 

138 

5 

3.68  - 2.23 

119 

6 

2.23  - 1.35 

121 

7 

1.35  - 0.82 

82 

8 

0.82  - 0.50 

76 

9 

0.50  - 0.18 

60 

10 

0.18  - 0.07 

10 

11-24 

En  < -1 

- 

